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Abstract 

We study the structure of the electronic states and the transport 
properties of a Kronig-Penney model with weak compositional and 
structural disorder. Using a perturbative approach we obtain an ana- 
lytical expression for the localisation length which is valid for disorder 
with arbitrary correlations. We show how to generate disorder with 
self- and cross-correlations and we analyse both the known delocalisa- 
tion effects of the long-range self-correlations and new effects produced 
by cross-correlations. We finally discuss how both kinds of correlations 
alter the transport properties in Kronig-Penney models of finite size. 



Pacs: 73.20. Jc, 73.20.Fz, 71.23.An 



1 Introduction 



The recent surge of interest in one-dimensional disordered models can be 
attributed in large part to the realisation that long-range correlations of 



the random potentials can significantly alter the structure of the electronic 
states, enhancing or suppressing their localisation in continuous intervals of 
energy [H [2j [3J SJ [5] . This discovery has paved the way for the construction of 
one-dimensional electronic, optical or electromagnetic devices with anoma- 
lous transport properties. Specifically, it is now possible to design filters with 
almost perfect transmission or reflection in pre-defined energy windows. This 
possibility, first foreseen at the theoretical level [21 E], was later confirmed 
experimentally [HE]. 

The periodic Kronig-Penney model was originally introduced in the early 
1930s to analyse the electronic states and the energy bands in crystalline 
structures [6]. In the 1980s the model was used to describe semiconduc- 
tor superlattices (see [7] and references therein). More recently, aperiodic 
variants of the Kronig-Penney model were used to analyse the transmission 
properties of a waveguide with long-range correlated compositional (U [5] or 
structural [8] disorder. The wide applicability of aperiodic Kronig-Penney 
models makes desirable to analyse in full detail the structure of their elec- 
tronic states. Two kinds of disorder are typically considered in the aperiodic 
Kronig-Penney model: structural and compositional disorder. In the first 
case the positions of potential barriers are randomised; thus, they do not 
coincide with the sites of the underlying lattice. In the second case the dis- 
order is introduced by assigning random heights or widths to the potential 
barriers. 

In this paper we study a Kronig-Penney model with delta-shaped po- 
tential barriers with weak correlated disorder of both the compositional and 
structural type. This model was first analysed in [9], where we presented 
an analytical expression for the inverse localisation length valid for any kind 
of self- and cross-correlations of the two disorders. Here we give a detailed 
derivation of this result, together with an extended discussion of various phys- 
ical implications. In particular, we describe how to construct two sequences 
of random variables, representing the positions and strengths of the barriers, 
with arbitrary self- and cross-correlations. We show that specific long-range 
self-correlations can create a delocalisation-localisation transition (within the 
limit of the second-order perturbative approach) and we analyse new effects 
produced by the cross- correlations of the two kinds of disorder. We then 
discuss how the both kinds of correlations affect the electronic transport 
properties in finite samples. 

This paper is organised as follows. In Sec. [2] we define the model under 
study and derive an analytical expression for the electronic localisation length 
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using the Hamiltonian map approach. In Sec. [3] we show how to construct 
sequences of correlated random variables for the positions and strengths of 
the barriers. We also demostrate how the localisation/delocalisation effects 
produced by long-range self-correlations of the disorder are modified by the 
cross-correlations. In Sec. H]we analyse how these effects alter the transport 
properties of finite-size models. The conclusions are drawn in Sec. 

2 Structure of the electronic states 
2.1 Definition of the model 

Here we define the Kronig-Penney model with weak compositional and struc- 
tural disorder. The model describes the motion of an electron in a succession 
of delta-shaped barriers; disorder is introduced by assuming that both the 
strength and the position of each barrier are random variables fluctuating 
around their average values. The corresponding Schrodinger equation is 



with the electron energy given by the square of the propagation wavevector 
within the wells 



In Eqs. (DP) and ([2]) we have used units such that h 2 /2m = 1; we will stick to 
this choice throughout this paper. 

In the model ([1]) the delta-barriers are located at the points d n = na + a n , 
where the random shifts a n with respect to the lattice sites na represent 
the structural (or positional) disorder. The compositional (or amplitude) 
disorder, on the other hand, enters via U n = U + u n , with the random 
variables u n corresponding to fluctuations of the barrier strengths around 
their mean value U. 

We define the statistical properties of the model in terms of the strength 
fluctuations u n and of the relative barrier displacements A n = a n+1 — a n 
(which turn out to be more relevant than the absolute displacements a n 
themselves). We restrict our attention to the weak-disorder case, defined by 
the conditions 



oo 



(1) 




(2) 




(3) 
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and 



(Ajy<l and (K)U<1. (4) 

The condition ([3]) implies that the compositional disorder is weak in the 
sense that the fluctuations of the barrier strength are small with respect 
to their mean value. The weakness of structural disorder is defined by the 
condition which requires that the relative displacements of the barriers 
should be small on the length scales set by 1/q and 1 / yU. 

In the weak-disorder case the statistical properties of the model can be 
analysed in terms of the averages and of the binary correlators of the random 
variables u n and A n . Both variables are assumed to have zero mean. In what 
follows we consider the normalised binary correlators 

(u n A n+k ) 



Xs(k) 



as given functions. We do not attribute any specific form to the correla- 
tors (JHJ); we only suppose that they are even functions of the index difference 
as a consequence of the assumed spatial isotropy and homogeneity in the 
mean of the model. 



2.2 The Hamiltonian map 



The electronic states of the Kronig-Penney model (CQ) can be analysed by 
means of the Hamiltonian map approach. This approach is based on the 
mathematical identity of the Schrodinger equation ([1]) with the dynamical 
equation of a classical oscillator with noisy frequency [TOl [11] . The stochas- 
tic oscillator corresponding to the Kronig-Penney model is defined by the 
Hamiltonian 



2 2*' 



U n 8(t — na + a r 



x 



(6) 



with the electron wavevector q playing the role of the unperturbed frequency. 
The mathematical analogy between the Schrodinger equation ([T]) and the 
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dynamical equation of the parametric oscillator (E]) allows one to analyse the 
behaviour of the electronic states of the Kronig-Penney model in terms of 
the trajectories of the kicked oscillator. This dynamical approach implies 
that the Schrodinger equation is solved as an initial-value problem; hence 
the Hamiltonian map method is equivalent to the transfer matrix technique. 

Given the Hamiltonian fl^]), one can integrate the corresponding dynami- 
cal equations over the time interval [d~, d~ +1 ) between two kicks. In this way 
one obtains the Hamiltonian map 



x n+l 
Pn+l 



•En 
Pn 



(7) 



with the transfer matrix 

cos [q (a + A n )] + (U + u n ) - sin [q (a + A n )} ~ sin [q (a + A n )] 

— q sin [q (a + A n )] + (U + u n ) cos [q (a + A n )] cos [q (a + A n )] 

(8) 

To analyse the trajectories of the Hamiltonian map (JTj), we follow the ap- 
proach proposed in [3J. Specifically we perform a canonical transformation 
(x n ,p n ) —> (X n ,P n ) such that the unperturbed motion reduces to a simple 
rotation in the new variables. The appropriate canonical transformation has 
the form 

X n ( M 

Pn 



Pn 



(9) 



with 



M 



a cos — 



V 



— qa sin ■ 



qa 



1 qa \ 

— sin — ' 

qa 2 
1 qa 

— cos — 
a 2 



where the parameter a is defined by the relation 

U 



a 



^ sin (qa) — — [cos (qa) — 1] 

sin (qa) [cos (qa) + 1] 

2q 



(10) 



Note that the canonical transformation (Q rescales time so that the new 
coordinate X n and the corresponding conjugate momentum P n have the same 
dimensions. 
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Since we are interested only in the weak-noise case, we can expand the 
new transfer matrix T n = M _1 T n M and keep only terms up to the second 
order. The Hamiltonian map then takes the form 



X n +i 

Pn+l 



= Tf + T W + T< 2 > + o(2) 
where, as desired, the unperturbed matrix represents a rotation 

V — sin (ka) cos (ka) 
with the rotation angle ka being defined by the equation 

cos (ka) = cos (ga) H sin (ga) . 

2q 



(11) 



(12) 



In terms of the Kronig- Penney model, k is the Bloch wavevector, while equa- 
tion (fl2"|) defines the band structure of the model. The explicit form of the 
first- and the second-order terms of the map (1111) is 



f(i) 



and 



X(2) 



/ 
\ 



sin (qa) 1 . , , , . 

— u n sin (ka) A r 

2q cr 



a- 



2 2 



1 + cos (qa) 



1 — cos (qa) . . 
2 u„ - cos (fca) A n 



« n — cos (ka) A r 



2q 
sin (qa) 



( cos (ga) + 1 _ 1 ?2 cog A 2 



2q 

sin (ga) 
2qa 



u n — a q sin (fca) A n 



, u n A n - ^g 2 sin (A;a) A 2 



/ 



ga 2 . , . . 1 2 . /, \ a 2 cos(ga) 
\ — =- sm (? a ) u n\ t + -g sm (fca) A„ ^~ 



t A n - -g 2 cos (ka) A 2 j 



A remark is in order here. The canonical transformation is well- 
defined for every value of the rotation angle ka except for the critical values 
ka = and ka = ±7r. In terms of the Kronig-Penney model, this means that 
our approach fails at the centre and at the edges of the first Brillouin zone. 
This can be seen by considering that Eq. ( !T2|) implies that 



sin (ka) 



sin (ga) — — [cos (ga) + 1] 
2g 



sin (ga) — — [cos (ga) — 1] 
2q 



This identity shows that, whenever the sine of the rotation angle is equal 
to zero, the a parameter ( imp either vanishes or diverges. In both cases the 
canonical transformation fl9]) is not properly defined, because some elements 
of the matrix M diverge. However, we would like to stress that the canon- 
ical transformation is perfectly well-defined as long as the rotation angle is 
arbitrarily close (but not exactly identical) to the singular values ka = and 
ka = ±7r. Hence our approach works well in every neighbourhood of these 
critical points. 

To analyse the evolution of the dynamical system ( ITTj) , it is convenient to 
switch from Cartesian to action-angle coordinates, defined via the equations 



X n = V2Xsin0 n 
P n = A/2Xcos6» n ' 

In terms of the new variables, the Hamiltonian map ( ITT1) takes the form 
T — n 2 T 

On+i = 9 n + ka - i [1 - cos (20 n + ka)] u n + ]- [C - cos (26 n + 2ka)} A r 
1 

+ - [2 sin (20„ + ka) - sin (40 n + 2ka)\ u\ 

~ ] 2 



+ - [2C sin (20 n + 2ka) - sin [A0 n + 4ka)] A„ 
+ — [sin (ka) — 2 sin (20 n + 2ka) + sin (A0 n + 3ka)] u n A n , 

with the ratio of the action variables being equal to 



(13) 



D 2 n = 1 + sin (20 n + ka) u n - sin (20 n + 2ka) A n 

+ ~[l-cos(20 n + ka)}u 2 n + -[l-(cos(20 n + 2ka)}A 2 n (14) 

+ [cos (20 n + 2ka) — cos (ka)] u n A n . 

In Eqs. ( fl3|) and (fl4|) we made use of the rescaled random variables 

sin (qa) , x U . 

= — : — — T u n and A n = , A n (15) 
g sm (/ca) sin (fca) 

and we introduced the short-hand notation 



gsin(/ca) 



C7 



<?a 2 + 



ga 2 



(16) 



We remark that the angle variable evolves independently of the action vari- 
able. 
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2.3 The localisation length 

This subsection is devoted to the task of evaluating the localisation length in 
the Kronig- Penney model ([TJ). The inverse localisation length is defined as 

N i 



n=l 



A. 



which, in dynamical terms, is equivalent to the Lyapunov exponent of the 
Hamiltonian map (ITT]) , i.e., 

A = lim — Y log ( ^] = —(\ogD 2 n ). (17) 

n=l v 7 



To compute the Lyapunov exponent ([17]) we expand the logarithm; within 
the second-order approximation one obtains 

A = — /sin (26 n + ka) u n — sin (26 n + 2ka) A n 
2a \ 

+ - [1 - 2 cos (26 n + ka) + cos (A6 n + 2ka)\ ul 

f ~, ( 18 ) 

+ -[1-2 C cos (26 n + 2ka) + cos (A6 n + Aka)} A 2 n 

— - [cos (ka) — 2 cos (26 n + 2ka) + cos (A6 n + 3ka)} u„ A n 

The terms in the angular brackets of the right-hand side (rhs) of Eq. ( ITS]) 
must be averaged over the angle variable with an appropriate distribution 
p(6). Typically, one uses a flat distribution p(6) = 1/(2tt). In fact, it is easy 
to see that in the absence of disorder the Hamiltonian map ( IT3]) reduces to 
the form 

Jn+l Jn 

6 n +i = 6 n + ka. 

Hence the angle variable has a fast dynamics compared to the action variable 
and quickly assumes a uniform distribution of values in the interval [0 : 27r] . 

This reasoning fails, however, when the rotation angle is either zero or 
a rational multiple of tt. In the latter case the noiseless angular map has 
periodic orbits which manifest themselves in the form of a slight modulation 
of the invariant measure when a weak noise is switched on. In particular, 
when one has 

ka — — 
n 
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the uniform distribution of the angular variable is modified by a perturbative 
term proportional to cos(2ra#) or sm(2n8). This modulation of the invari- 
ant measure alters the value of the inverse localisation length flTHj) only for 
n = ±1, and n = ±2, because only second- and fourth-order harmonics are 
present in the rhs of Eq. f|T8|) . The values ka = and ka = ±7r correspond 
to the centre and the edges of the first Brillouin zone, i.e., to the edges of the 
energy bands. The values ka = ±7r/2 of the Bloch vector, on the other hand, 
correspond to values of the energy close to the band centre. For these special 
values of the rotation angle the assumption of a flat invariant distribution 
must be abandoned and the specific form of p{6) has to be determined before 
one can compute the average in Eq. (I18p . The non-uniform distribution of 
the angular variable produces anomalies in the localisation length; we will 
analyse them elsewhere. 

Except for the special cases mentioned above, one can compute the av- 
erages in Eq. f|T8|) with a uniform angular distribution. The corresponding 
result is 



A 



1 



+ 



1 
2a 



a l 



vf n ) + (A„) - 2(u n A n ) cos (ka 
u n sin (28 n + ka)) — (A n sin (2Q n + 2ka)) 



(19) 



To proceed further, we need to compute the noise-angle correlators which 
appear in formula (jT9l) . One can generalise the method used in [2] and 
introduce the correlator between the kick strength and the angle 

n = (w n exp (i26 n _i)) 

as well as the correlator between the kick timing and the angle 

si = (A„exp (i29 n ^i)). 

Both correlators satisfy recursive relations, which can be derived as follows. 
Keeping only second-order terms in the disorder strength, one can write r/_i 
in the form 

r/_i = (u n exp (i29 n -i + i) J = (u n exp (z20 n _/) exp (i2ka) 

x 1 1 - i [1 - cos (26>„_; + ka)] u n -i + cos (20^-z + 2ka)] A n _; 

(20) 
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For weak disorder the triple and quadruple noise-angle correlators can be 
factorised and averages over the angle variable can be carried out with a 
uniform distribution. The result is 

ri-i = rjexp (i2ka) + -exp (ika) (u n u n -i) - -(u„An-i)- 

Multiplying both members of this identity by exp [i2ka (I — 1)] and summing 
over / from zero to infinity, one obtains 



r = - exp (—ika) ^^(u n u n _i) exp (i2kal) — - exp (—i2ka) V^(u n A n _j) exp (%2kc 
1=1 1=1 



From this expression one can easily compute the correlator between the angle 
and the kick strength in Eq. (|T9|) ; in fact, one has 

(it„sm (26 n + ka)) = Im [r exp (ika)} 

1 °° 1 °° 

= -^(u n u n -i) cos (2kal) - -^(u n A n _i) cos [ka (21 - 1)]. ^ > 



i=i 



i=i 



Following the same approach, one obtains that the correlator between the 
angle and the timing of the kick is 

(A n sin (26 n + 2ka)) = Im [sq exp (i2ka)} 
1 °° 1 °° 

i=i i=i 

Plugging the noise-angle correlators f l2"Tj) and (|22p in Eq. ( 1T§|) . one obtains 
that the Lyapunov exponent is 

A = UuDW^ka) + (A 2 n )W 2 (ka) - 2(u n A n ) cos (ka) W 3 (ka)} (23) 

where the functions Wi(ka), defined by the identities 

(u n u n+ i) 



W x (ka) 
W 2 (ka) 
W 3 (ka) 



1 + 2^ 



i=i 

OO 



■ cos(2kal) 



1 , oV^\^ n ^ n +'/ rm a 
1 + 2 > — — — — cos(2kal) 

i=i ' n > 

1 i oV^ \ U nA n +l) / n7 A 

1 + 2 > — — ^— p cos(2kal), 



(24) 
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{u n A n ) 



are the Fourier transforms of the normalised binary correlators §J5§ . Eq. f)23|) , 
first derived in [9], shows that the Lyapunov exponent is the sum of three 
terms, with the first two addends describing the effects of purely composi- 
tional and structural disorder, while the third term is due to the interplay 
between these two kinds of randomness. 

We stress that expression (1251) is valid for every value of the energy within 
the allowed energy bands, with the exception of small neighbourhoods at the 
band edges (i.e., for ka ~ and ka ~ ±7r) and close to the band centre (i.e., 
for k ~ -k /2a). 

3 Designed mobility edges 

3.1 Generation of self- and cross-correlated disorders 

In this section we discuss how specific long-range self-correlations of the dis- 
order can weaken or enhance the localisation of the electronic states and 
produce pre-defined mobility edges. We then analyse the additional effects 
of the cross-correlations between compositional and structural disorders. For- 
mula ( 1231) shows that the localisation length diverges in the energy intervals 
where the power spectra vanish. The problem, therefore, is to generate 
two successions of random variables {u n } and {A n } with self- and cross- 
correlators such that the power spectra (I24|) vanish in pre-assigned energy 
windows. We are dealing with an "inverse" problem which, as such, has no 
unique solution. Here we propose one, defining an algorithm for the con- 
struction of two successions {«„} and {A n } with the required features. 

As a first step, we consider two successions {X^} and {X { n ] } of inde- 
pendent random variables with zero mean and unit variance. In other words, 
we require that 



for i,j = 1,2 and n, m G Z. In terms of these variables, we can construct 
two inter-correlated successions 



with r] being a real parameter which determines the degree of inter-correlation 



(Xg>) = Q and (Xg>X%>) = S ij 6 i 




Xn cos r] + Xn ^ sin 7] 

Xn Sin 7] + Xn^ COS T) 



(25) 
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of the Y variables. In fact, one has 

= and (y«yW) = 5nm [ 6ij + (1 - <y sin (2^)] 

and this shows that the natural range of variation of the parameter 77 is the 
interval [— 7r/4, 7r/4] with r] = tt/A corresponding to total correlation of the 
in and F n variables and 77 = —tt/A to total anticorrelation (for rj = no 
cross-correlations exist). 

The next step consists in "filtering" the cross-correlated white-noise suc- 
cessions and {Fi 2 "*} in order to obtain two cross- and self-correlated 
sequences {«„} and {A n }. For this purpose we express the u n and A n random 
variables as convolution products of the form 

00 00 
u n = a * Y n-t and A n = ^ Y n%- (26) 

k=— oo k=— oo 

and we look for appropriate successions of the coefficients ot k and f3 k - We 
suppose that ot k = a.- k and (3 k = 0-k- Starting from Eq. ( )2"6"|) one obtains 
that the average values of the variables u n and A n vanish 

(u n ) = and (A n ) = 

and that the binary correlators take the values 

00 

(u n u n+k ) = ^2 aiai+k, (27) 
/=— 00 

00 

(A n A n+fe ) = AA+fc, (28) 

i=— 00 

00 

(A n u n+k ) = (u n A n+k ) = 22 a iPi+k sin (2r]) . (29) 

i=— 00 

Comparing Eqs. (I27p and ( 128|) with the binary self-correlators in Eq. (jSJ), one 
arrives at the equations 

Y a i a i+k = (u 2 n )xi{k) 

l=—oo 

00 

X>A+* = (An)X2(fc) 
!=— 00 
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whose solution is 



2 rl 2 

«fc = — / y/ {v%)Wi(x) cos (2kx) dx, 

2 J ° r /2 (30) 
(3 k = - y/(Al)W 2 (x) cos (2fcp) dx. 
^ Jo 

Inserting the coefficients (1301) in the convolution products fl26|) one obtains 
random variables {u„} and {A n } with zero average and pre-assigned self- 
correlators Xi(^) an d X2(^), as required. To complete the picture, one 
should specify the inter-correlation parameter rj which, together with the 
coefficients f[3T)j) defines the cross-correlator (J2"5j) . 

Thus, we have obtained a recipe that gives two random successions {u n } 
and {A n } in terms of the pre-defined normalised self-correlators Xi(^) an d 
X2(k), or, equivalently, of the pre-assigned power spectra Wi(ka) and W2(ka) 
and of the inter-correlation parameter rj. We can now see what is the form 
of the inverse localisation length (1231) when the random variables u n and A n 
are generated with the above- described method. After some algebra, one can 
write the cross power spectrum as 

(u n A n ) 

Inserting this result in expression (123]) one obtains 
1 



(u 2 n )Wi(ka) + (A 2 n }W 2 (ka) 



2\J(ul) (Al)Wi(ka)W 2 (ka) cos (ka) sin (2?y) 



(31) 



3.2 Numerical examples 

To test the validity of formula ( 13~T|) we construct a pair of random sequences 
{u n } and {A„} which enhance the localisation of the states in two energy 
windows [qf, q%\ and g|] and delocalise the states in the rest of the band. 
For the sake of simplicity we consider compositional and structural disorders 
with identical self-correlators of the form 

Xi{n) = Xa(n) = —r. ^ t sin i 2k 2an) - sin (2/^cm)] (32) 

2 (k 2 — k\)an 
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The power-law decay of this function is the signature of long-range auto- 
correlations of the disorder. The correlators (1321) correspond to the power 
spectra 

— ■ if k G \ki,k 2 ] 

2a(k 2 -k l ) [ 21 . (33) 

if he [0, h] U [k 2 , |J 

In Eqs. (1321) and ( 1331) . k\ and k 2 represent the two Bloch vectors corresponding 
to the mobility edges q\ and q 2 , i.e., 

q(ki) = q 1 and q(k 2 ) = q 2 - 

Note that the power spectra (1241) are periodic functions of period n/2a and 
therefore it is enough to define them in the interval [0, 7r/2a]. 

For our first numerical example we chose disorder strengths equal to 
a/ (u^) = a/ (Af) = 0.04 and a mean field value U = 0.7. We considered 
two Bloch vectors k x = n/5a and k 2 = 2n/5a, which, for the given value of 
U, correspond to mobility edges at 

qi = 0.327, q 2 = 0.476, q 3 = 0.652, q A = 0.838 

with ^ = qi/ir. Note that, with our choice of U, the propagation wavevector 
q within the first energy bands spans the interval ranging from qi/ir = 0.259 
to q r /n = 1.0. Fig. [TJ shows that the theoretical predictions of Eq. ( 1311) 
are well-matched by the numerical results. As expected, the long-range self- 
correlations ( 132]) produce sharp mobility edges, whereas the cross-correlations 
can substantially change the spatial extension of the localised states. In 
Fig. ([T]) the case without cross- correlations (rj = 0) is compared with the 
two extreme cases of total positive (rj = vr/4) and negative (rj = — 7r/4) 
cross-correlations. 

As a second example, in Fig. [2] we represent the theoretical predictions 
and numerical results for the complementary case, characterised by the same 
mobility edges of the previous example but with inverted windows of localised 
and delocalised states. (The values of the disorder strengths and of the 
mean field are the same as in the previous case.) We observe again a good 
correspondence between the numerical values and the inverse localisation 
length (|3Tj). 

To conclude this section, we observe that the spatial extension of the 
localised states can be significantly reduced by squeezing the energy windows 
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Figure 1: Inverse localisation length versus q/ir. Lines correspond to the 
predictions of Eq. ( 13TT) ; the symbols to numerical results. 



of localised states. This localisation enhancement, first discussed in [5], is a 
consequence of the normalisation condition 

t/(2o) 7r 
Wi(ka)dk = — , 

which follows from the fact that Xi(0) = 1. To illustrate this effect, we 
consider a third example, with all parameters unchanged with respect to 
the first one (U = 0.7, a/T^U = \l (^n) = 0.04) except for the size of the 
windows of localised states, which we shrink by setting k\ = 0.29n/a and 
hi = 0.3l7r/a. In terms of the propagation wavevector, these values of the 
Bloch vectors correspond to mobility edges at 

q x = 0.388, q 2 = 0.402, q 3 = 0.736, g 4 = 0.755. (34) 

The theoretical and numerical results for the inverse localisation length, rep- 
resented in Fig. [3], display a clear enhancement of localisation in the narrowed 
energy intervals of localised states. The stronger localisation magnifies the 
effect of the cross-correlations of the compositional and structural disorder, 
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Figure 2: Inverse localisation length versus q/it. Lines correspond to theo- 
retical predictions and symbols to numerical results. 



enlarging the difference between the extreme cases of total positive and neg- 
ative inter-correlations. 

4 Transport properties 

In the previous sections we have considered the structure of the electronic 
states in a infinite Kronig-Penney model. We now turn our attention to 
the relevant problem of electronic transmission through finite disordered seg- 
ments. We consider the case of a random Kronig-Penney model extending 
over N lattice sites sandwiched between two semi-infinite perfect leads. From 
the mathematical point of view, this means that the variables u n and A n in 
the Schrodinger equation (TjQ) are defined as before for 1 < n < N, but vanish 
for n < 1 and n > N . We assume that the left lead carries an incoming and 
a reflected wave with wavevector q, 

^ = e iqan + re -iqan for n = 0, —1, —2, . . . 
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Figure 3: Inverse localisation length versus q/ir. Lines correspond to theo- 
retical predictions and symbols to numerical results. 



while a transmitted wave propagates in the right lead 

^ n = te iqan for n = N +1,N + 2,... 

To determine the transmission coefficient T(q) = \t(q)\ 2 we follow the 
variant of the transfer-matrix approach introduced in [12j . After eliminating 
the momenta from the map (j7|), with the obvious substitution x n — > ip n we 
obtain the equation 

f^n+l + , r . 1 rrV'n-l 



sin [q (a + A n )] sin [q (a + A n _i)] ^ 

cot \q (a + A n )] + cot [q (a + A n _i)] + i (U + u n ) | i> n , 

which defines the tight-binding system corresponding to the Kronig-Penney 
model ([1]). Note that the compositional disorder in the Kronig-Penney model 
manifests itself as diagonal disorder in Eq. fl35|) . whereas the structural dis- 
order shows up in both the diagonal and off-diagonal random coefficients of 
the tight-binding model (|35p . Using the transfer-matrix formalism, one can 
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write Eq. (I55|) in the form 



^ n+1 ^ - M n ( f n ) (36) 



with the elements of the matrix M n being 

(M„) n = |cot[g(a + A n )]+cot[g(a + A n _ 1 )]+i(C/ + « n )}sin[ g (a + A n )] 
riu \ _ sin [q (a + A n )] 
sm [g (a + A n _i)J 

(Mn) 21 = 1, 

(M n ) 22 = 0. 

The total transfer matrix connecting the wavefunctions in the two leads is 
then 

M = M JV M JV _i---Mi. (37) 

Following Pichard [13] , one can write the transmission coefficient in terms of 
the elements of the matrix (I3TI) . One thus obtains the expression 

T N{q) = 4sm2 ^( detM ^ - 2 . (38) 

|M 21 - Mi 2 + Mjcje* - Mne-^l 2 

We remark that the factor (det M) 2 enters the previous formula because the 
determinant of the transfer matrices M n is not unitary. More precisely, one 

detM= Sin[ r 9( ; + ^ ] . 

sm [q (a + A )J 

To avoid the nuisance of a total transfer matrix with non-unitary determi- 
nant, one can redefine the random barrier as including an extra site with no 
disorder on each side. From now on we will adopt this convention and drop 
the determinant factor in formulae derived from Eq. fl38|) . 

To express the transmission coefficient ( 1381) in a more convenient way, we 
perform a similarity transformation 

Q n = RM n R 1 

with 

R= ( 1 ° 

cos (qa) j sin (qa) — l/sin(ga) 
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and put the map in the equivalent form 



, On n ■ (39) 



This allows one to cast the transmission coefficient fl38|) in the more appealing 
form [12] 

= ^- A — , (40) 

2 + vi • vi + v 2 • v 2 

with Vi and v 2 being the two-component vectors defined as 

Vl = Qtv+iQtv • ■ ■ QiQo 

and 



v 2 — Qv+lQjV ' - - QlQc 



Eq. ( l40l) is very convenient for the numerical evaluation of the transmission 
coefficient, which can be efficiently determined by computing the evolution 
of the initial vectors (1, 0) and (0, 1) under the map ( [391 . 

Using this approach, we numerically investigated the effects of the long- 
range correlations of the disorder on the transmission coefficient T^(q). In 
the absence of cross-correlations, only two length scales are relevant for the 
transport properties of the disordered segment, namely, the length of the 
segment itself, L = Na and the spatial extension Zi oc of the localised states 
in the infinite random model. In principle one should also consider the lo- 
calisation length /del of the delocalised states, because the vanishing of the 
second-order Lyapunov exponent fl23l) does not necessarily imply that the 
higher-order terms must also be zero. However, for weak disorder the spatial 
extension of the delocalised states is so large that it can be considered as 
infinite with respect to the length L of the random barrier. 

In the localised regime, i.e., when the condition /i oc <C L is fulfilled, our 
numerical experiments show that the windows of delocalised states created 
by long-ranged self-correlations of the disorder survive in the case of dis- 
ordered segments, where they manifest themselves as energy windows with 
transmittivity close to one. This agrees with the results already obtained, 
both numerically and experimentally, for barriers with purely compositional 
or structural disorder (3J |8]. Our numerical data also show that the cross- 
correlations of the disorder have additional effects on the transport proper- 
ties. Their influence is most evident in the localisation windows created by 
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long-ranged self-correlations of the disorder, where cross- correlations either 
enhance or diminish the electronic transmission according to whether they 
increase or decrease the localisation length. 

The effects of both the self- and the cross-correlations of the disorder can 
be appreciated in Figs. H] and |5l which represent the transmission amplitude 
for two lengths of a random sample with the same disorder characteristics 
of the third example analysed in Sec. 13.21 In other words, we consider a 
random Kronig- Penney barrier with mean field U = 0.7, disorder strengths 
■sj (11%) = \J (Af) = 0.04, and windows of localised states [91,52] and [93,94], 
with edges defined by Eq. (I34p . Each figure represents the transmission 
coefficients for the cases of total positive and negative cross-correlations. 



Positive cross-correlation 
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Figure 4: Transmission coefficient T N (q) versus q/n for random sample of 
N = 100 sites. 

When the cross-correlations are taken into account, it is better to char- 
acterise the spatial extension of the localised states in terms not of one but 
of two length scales, i.e., /j£ c (?) and l^ c (q), which correspond to the maxi- 
mum and minimum values of the localisation length obtained by varying the 
inter-correlation of the compositional and structural disorders. In the present 
case, one has l^ c ~ 100a and ~ 300a in the localisation window of lower 
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Negative cross-correlation 




Figure 5: Transmission coefficient T/v(g) versus q/% for random sample of 
N = 200 sites. 



energy and ~ 200a and ~ 400a in the higher-energy localisation win- 
dow. The data represented in Figs. H] and [5] were obtained for the sample 
lengths L = 100a and L = 200a, which correspond to the intermediate region 
between the ballistic regime, characterised by L <C and the exponentially 
localised regime for which Zj£ c ^ L. 

Two main features emerge from the numerical data. On the one hand, as 
L is increased the mobility edges become sharper and sharper, because the 
localisation of the electronic state in the windows [gi, g 2 ] and [g 3 , q^} becomes 
more and more effective. On the other hand, the existence of two localisation 
lengths implies that the regime of exponential localisation is attained with 
different speed according to the type and degree of cross-correlations. This 
is the physical origin of the difference in the transmission coefficients for 
positive and negative cross-correlations. 

In general terms, one should notice that, when the length of the sam- 
ple is less than the minimum localisation length / 1 ~ c , the localisation of the 
electronic states is very weak, and this can produce strong sample-to-sample 
fluctuations of the transmission coefficient which partially mask the effects 
of cross-correlations. As the length of the barrier is increased, these effects 
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become clearly discernible and their influence peaks in the intermediate re- 
gion between the ballistic and the localised regimes. If one lets the length 
of the random sample grow further, the regime of exponential localisation is 
eventually reached for both kinds of cross-correlations. In this regime both 
transmission coefficients are exponentially small in the localisation windows 
and their absolute difference is therefore hard to detect. 

5 Conclusions 

In this paper we have analysed the structure of the electronic states and 
transport properties of a Kronig-Penney model with weak compositional and 
structural disorder. We have shown that specific long-range self-correlations 
of the disorder can enhance or suppress the localisation of the electronic 
states in pre-defined energy windows and that an additional modulation of 
the localisation length can be obtained by cross-correlating the two kinds of 
disorder. A method to generate compositional and structural disorder with 
arbitrary self- and cross-correlations has been given. 

The analysis performed in the second part of this work shows that, by 
carefully selecting the statistical properties of the disorder, one can produce 
random samples with peculiar transport properties. In particular, it is pos- 
sible to obtain almost perfect transmission in pre-defined energy windows. 
Although self-correlations play a dominant role in shaping the transport prop- 
erties of finite-size random chains, cross-correlations can be used to fine-tune 
the transmission coefficient. 
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